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ABSTRACT 

Pre-flood phytosoeiologieal measurements were made in eight apparently mature riparian forest 
eommunities along the upper Guadalupe River in eentral Texas. After a high magnitude 100+ year flood 
in 1978, six eommunities were revisited and post-flood phytosoeiologieal measurements were made to 
examine flood effeets on the riparian forest eommunities. Some riparian forests were eompletely 
destroyed with woody plant trunks snapped or plants uprooted and washed downstream. We viewed 
these as early sueeessional eommunities, and did not study them. The flood eaused ehanges in the 
riparian forest eommunities we studied previously ineluding reduetions in eommunity basal area and 
density, as well as speeies basal area, density, and pereent eover. Overall, eommunity speeies 
eomposition was unehanged based on Shannon-Weaver diversity values (H’). Speeies with the greatest 
reduetion in eommunity basal area or density were Taxodium distichum (Bald eypress), Carya 
illinoinensis (Peean), Ulmus americana (Ameriean elm), Celtis laevigata (Haekberry) and Acer negundo 
(Box-elder). Greatest eommunity damage seemed to be related to tree size and bank position. 
Apparently, the larger the tree or speeies, the greater the damage sustained in the flood. Larger trees or 
speeies at the edge of the river sustained the greatest eommunity ehange while the smallest ones at 
greater distanee from the river edge had the least ehange. However, aerial photographie evidenee of 
major river bends, areas that we did not study, suggested that water was so deep and the flood foree so 
great that it eut aeross the bend removing most trees and vegetation in its path. Published on-line 
WWW.phytologia.org Phytologia 99(4): 226-237 (Dec 18, 2017). ISSN 030319430. 
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Understanding eommunity strueture and faetors that determine strueture is a diffieult task. 
Usually more than one foree or faetor is involved and they eould be biotie, abiotie or a eombination. 
Abiotie forees or physieal forees sueh as fires, wind, tidal waves, voleanoes, earthquakes, floods and 
landslides ean have major impaets on plant eommunities and their eomponent speeies (Begon et al. 
2006). Biotie forees or faetors sueh as herbivory or eompetition ean have similar eommunity effeets and 
usually more than a single foree is in play (Begon et al. 2006). All of these forees or eombinations ean be 
important in direeting eommunity sueeession (Whitmore 1974; Putz et al. 1983; Begon et al. 2006; Van 
Auken and Bush 2013). 

Riparian forest eomposition and strueture has been studied for many years, ineluding studies in 
eentral Texas (Van Auken et al. 1979; Ford and Van Auken 1982). There are spatial gradients that at least 
partially eontrol the speeies eomposition of riparian forests. For example, flood tolerant trees or those 
requiring more water are elose to the river or stream bank and those less tolerant are more distant from the 
water edge (Bush and Van Auken 1984). Effeets of prolonged inundation on tree growth for speeies on 
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low gradient streams has been well investigated (reviewed by Hupp 1988; also see Johnson et al. 2000). 
Short tenn effects of deep, rapidly moving water on woody species and community structure are not as 
well understood because they occur so mlfequently. Some qualitative studies of flood damage have 
appeared m the literature mcludmg mention of complete destruction of all trees in parts of floodplain 
forests, as well as uprooting, breaking and shearing of stems (Shull 1944; Gardner 1951; Eisenlohr 1952; 
Hack and Goodlett 1960; Lindsey et al. 1961; Sigafoos 1964; Hallgren 1979; Sullivan 1983). 

Recent studies in central and western North America have demonstrated specific flood effects 
including tlie importance of floating debris, logs and whole trees (Johnson et al. 2000). In addition, 
changes m riparian forests on ephemeral streams have been examined (Friedman and Lee 2002). 
Predictive theories have also been developed that consider various aspects of extreme flood events 
including water force, material transport, flood frequency and conditions of the floodplain or bank (Dixon 
et al. 2015). 

Early studies of riparian communities did not have much access to or use of aerial photography 
and overviews of damage to large areas were not easily accomplished, but that has changed (Figure 1, 
Johnson et al. 2000; Friedman and Lee 2002; Dixon et al. 2015). These studies mamly examined flood 
effects on eaiJy establishmg or early successional floodplain species and communities. It has been 
suggested that temperate floodplain forests are good examples of communities that show a need for 
chronic, patchy disturbance (Horn 1976). The area of some of these patches may be extensive and show 
the results of regular destruction which would appear to be cyclic and triggered by periodic high 
magnitude floods. 

Southwestern and central Texas has one of the most varied rainfall-mnoff regimes in the United 
States (Patton and Balcer 1976). The rivers occur in bedrock canyons that are deeply entrenched in the 
Edwai'ds Plateau. The hills are steep, the limestone-derived soils are shallow and the vegetation is 
comparatively spai'se (Van Auken et al. 1979; Kochel 1988). Flood stages have been recorded between 
25 and 30 m above normal water levels in some areas (Figure 2, Kochel 1988). On August 1-4, 1978, 
more than 120 cm of rain from a stalled wealc tropical depression (Tropical Stonn Amelia) fell on the 
watersheds of tlie Medina and Guadalupe Rivers in central Texas (Hallgren 1979; Sullivan 1983). Most 
of this ram fell during a period of 24 hours causing a major flood. There were 27 deaths in central Texas 
attributed to the rains and subsequent floods. Normal rainfall for this area is about 81 cm/year, but is 
highly variable (Arbingast et al.l976; NOAA 2004). The return time for such a flood has been estimated 
to be several hundred years (Baker 1975; Hallgren 1979; Kochel 1988). Effects of the runoff on the 
riparian forest communities m some areas were extensive (Figure 1 and 2).Trees greater than 2 in in 
diameter were up-rooted or the stems were completely shattered with tlie stem being removed. Others 1 
m in diameter were found completely bent-over, parallel to the ground, with vertical splits in the stems 
(Hallgren 1979; Sullivan 1983). 


PURPOSE 

The purpose of this study was to identify and quantify changes in the woody vegetation of the 
upper Guadalupe River riparian forest caused by a high magnitude flood. An additional purpose was to 
show that the flood caused openings in the tree canopy reduced overstory community cover. The same 
pre-flood and post-flood riparian communities will be examined to demonstrate the vegetation changes. 

METHODS 

After analyzing and comparing results of an earlier phytosociological study of the Guadalupe 
River riparian forest (Ford and Van Auken 1982), we felt there was considerable variability in our data 
and additional effort would be necessary to understand these forest communities. Consequently, we 
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located and sampled eight additional inner floodplain forest eommunities to inerease our knowledge of 
this system. These riparian communities were relatively narrow (15-25 m) and relatively steep (slope ~ 
21°) usually with agrieulture land above the riparian woodland or forest. 

Pre-flood sampling. Durmg June and July 1978, phytosoeiologieal data was colleeted from 
eight riparian communities in Comal, Kendal and Kerr counties m central Texas. All study locations 
were above Canyon Reseiwoir, between approximately 50 and 75 mi north or northwest of San 
Antonio, Texas. We examined the distribution of species and community characteristics on the bank 
of the imier floodplain, between the river and the grassland or pasture using adjacent 3 x 100 m 
rectangular quadrats (Van Auken et al. 2005). Basal area, density, cover and species diversity 
(Shamion and Weaver 1949) were detennined and mean values were ealculated. Relative values for 
basal area and density were calculated but are not presented. Relative canopy cover was estunated by 
visual sightmg along a line perpendicular to a level surface at 25 m mtervals along the upper edge of 
each quadrat. Plant species names and identifications were detennined usmg “Vascular plants of 
Texas” (Conell and Johnston 1979) and eonfirmed at the United States Department of Agrieulture 
web site (USDA 2016). 

Post-flood sampling. Durmg September 1978, six of the original communities were re-sampled 
(after the flood). Because of community destruction and difficulties getting tlu'ough the flood debris 
(Figure 2), a different but equivalent sampling procedure was used (Cottam and Curtis 1956). The point- 
centered-quarter (PCQ) procedure was used to sample each community. Twenty-five points were 
sampled at 10 m intervals along a transect line parallel to and approximately 5 m from the river edge. A 
perpendicular to the transeet line was established creating four eomers where the perpendieular line 
erossed the transect line. The point was at the intersection of the two lines. For the woody plant closest 
to the point in each comer, the speeies was identified; then, point-to-plant distance, and eireumference at 
approximately 15 cm above ground level were measured and recorded. Only live plants were measured. 
Basal area, density, cover and species diversity were determined and mean values ealeulated. Relative 
values for basal area and density were calculated but not reported. Trees with sapped or broken stems 
were not counted or measured. 

Data analysis. Significant differences between the pre-flood and post-flood eommunities were 
determined. Total community mean basal area, density, cover and diversity were compared. Statistieal 
signifieance of eomparisons was detennined at the P=0.05 level using paired Mests (Sail et al. 2012). 

RESULTS 

Mean total community basal area was reduced 22% (Table 1). Four of the largest species, 
Taxodium distichum (L.) Rich. (Bald cypress), Cary^a illinoinensis (Wang., K. Koch., Pecan), Ulmus 
amehcana L. (American elm), and Celtis laevigata Willd. (Hackbeny) had average decreases in basal 
area of 8-33% (Table 1). The next four speeies in the table also had large losses in eommunity basal area 
(Table 1). The remaining speeies in table 1, low eommunity basal area speeies, had small or no 
reductions m basal area. 

Riparian forest community density was reduced 26% (Table 2). Caiya illinoinensis, the species 
with the highest pre-flood density, changed very little after the flood. Species that were found higher up 
on the floodplain had little or no change in mean density post-flood (C. illinoinensis and U. crassifolia). 
Species found at the river edge had the greatest reduetion in density (22-30%, Table 2). However, less 
eommon or small floodplain speeies had little or no ehange in mean density. 
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Species richness in the riparian communities was 31 woody species (Tables 2 and 3). Most of the 
species had low density and basal area. The three ti'ee species with the largest mean species basal area in 
descending order were T. distichum. C. illmoinensis and U. atnericana (Table 4). The largest tree 
measured was a Taxodium distichurn that had a diameter of 1.15 m. Relative cover was reduced 28% by 
the flood (from 96±4% to 69±5%, paired ^test, P< 0.05). Species diversity (H) did not change post¬ 
flood (2.06±0.34 pre-flood and 2.01±0.35 post-flood, paired /-test, R>0.05). 

DISCUSSION 

Biotic and abiotic forces can have major effects on plant communities (Whitmore 1974; Putz et 
al. 1983; Begon et al. 2006; Van Auken and Bush 2013). Changes in worldwide plant communities 
including North American and central Texas plant communities have been going on since prehistoric 
times. Changes in riparian communities have also occurred in the past and will continue to change in the 
future (Barbour and Billings 1988; Crawley 1997; Begon et al. 2006; Van Auken and Bush 2013), A 
successional sequence in these various riparian communities is involved, but fire, wind and herbivoiy/ do 
not seem to be the forces directing the succession. In the riparian woodlands, the major forces influencing 
the composition and structure of the communities is flood frequency, the time between events, the 
magnitude of the flood, river bank position and amount of floating debris (Johnson et al. 2000; Friedman 
and Lee 2002; Dixon et al. 2015). 

Curiously, we thought the riparian communities we studied pre-flood were mature coimmmities 
based on the literature available, the species present and our understanding of these communities at that 
time (Ford and Van Auken 1982). The riparian coimnunities we examined pre-flood did not include any 
Salix trees (willows) or Populus species (cottonwoods) and only a few Platanus occidentalis (sycamore). 
These species are all laiown or reported to be early successional riparian woodland species (Barbour and 
Billmgs 1988; Crawley 1997; Van Auken and Bush 1988). We have not found any studies concerning 
succession in the central Texas ripaiian forest communities. However, we speculated that T. distichurn is 
an early successional species (Ford and Van Auken 1982). 

Since that time, we have not found juveniles in any of the riparian communities that we 

surveyed (Van Auken 1988, 1993). Some preliminary studies suggested that T. distichurn is a high light 
requiring species (not found below the canopy), indicatmg that it is an early successional species (Van 
Auken and Taylor, unpublished). Another species found, Celtis laevigata, is a shade tolerant and high soil 
nitrogen requiring species (typically found below a canopy). Characteristics suggest it is a later 
successional species, one that would grow in the shade of an overstory canopy (Van Auken and Bush 
1985; Bush and Van Auken 1986; Van Auken and Bush 2013) and this cun'ent study suggests tliat as 
well. 

In central Texas, populations of upland species like live oak and Texas persimmon show 
evidence of recruitment (Van Auken et al. 1980). Conversely, other upland deciduous species like 
Quercus sinuate (Texas red oak), Q. glaucoides (Lacy oak), Pninis serotina (black cherry), Fraxinus 
texensis (Texas ash), and Juglans nigra (black walnut) had few or no new juveniles in the smallest size 
classes examined (Van Auken 1988, 1993). Therefore, there is little or no recruitment of juveniles of 
these species into the adult populations. 

This lack of recruitment of adults occurs in both upland cormnunities and riparian communities. 
Possible reasons for this lack of recruitment are changes in local climate, changes in the disturbance 
regime, increases in herbivore populations, few large carnivores, or episodic establishment that reflect 
local environmental conditions (Harper 1977; Harcombe and Marks 1978; Beschta and Ripple 2009; Van 
Auken and Bush 2013). Records do not indicate local clunate changes, although they are occunmg and 
will be more obvious in the future (Grunstra and Van Auken 2015). Proximate conditions may be more 
important such as populations of white-tailed deer in central Texas that have increased dramatically in the 
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last 75+ years (Doughty 1983). These large deer populations have been implieated in the laek of 
reeruitment of various upland deeiduous woody speeies in eentral Texas (Russell and Fowler 1999, 2002, 
2004; Nelson-Diekerson and Van Auken 2016). 

This lack of recruitment has been demonstrated for populations of T. distichum (Van Auken 
1988), which is similar to trends from floodplain forests in the eastern United States (Sharitz et al. 1990). 
Reasons seem to be high liglit levels that are not present below a canopy but would occur in disturbances 
(Neufeld 1983; Titus 1990). If T. distichum is to remain a component in these central Texas riparian 
forests, disturbances of these ecosystems should continue. Populations of C. laevigata m tliese same 
riparian communities seem to have large numbers of seedlings and there seems to be recruitment into the 
adult population. Celtis laevigata is shade tolerant and a mature riparian forest dommant species unlike 
T. distichum that appears to be shade intolerant (Van Auken and Bush 1985). 

Elevated levels of CO 2 and concomitant increases in temperature will undoubtedly effect the 
species in this area and the community structure in the future, but how is uncertain. Elevated levels of 
CO 2 have caused mcreased levels of photosynthesis in some species, but associated increases in 
temperature had no effect compared to cuiTent temperatures (Grunstra and Van Auken 2015). 

Riparian communities in central Texas, like other plant communities will change with 
disturbances, old communities with large/old T. distichum will mostly be destroyed by periodic high 
magnitude floods (depending on river position and water depth). Thus, the age of old T. distichum trees 
or communities in this area could be used to predict the approximate time for the next high magnitude 
flood event or perhaps when the previous one occurred, but neither has been done. There are probably 
other factors to be considered when maldng predictions, such as sand or gravel bars, chutes, and the cut 
banks of curves in the river that may be extremely susceptible to flood damage (Figure 1). We found no 
evidence of damage to the Guadalupe River riparian forest trees due to inundation as in other river 
systems, perhaps because floods of central Texas Rivers seem to occur and dissipate veiy' quickly (Baker 
1975; Hallgren 1979; Sullivan 1983). 

Larger trees seemed to experience the greatest damage, which suggests that they are the most 
susceptible to the force of the flood. The smallest species and individuals would be subject to the lowest 
water velocity, because of their size and position in the water column. Woody plants with small canopies 
would also offer the least resistance to the flood water. Debris trapping (Sigafoos 1964; Jolmson et al 
2000) would also be reduced because small individuals would be completely submerged and removed 
from the surface where the floating debris occurs. Larger trees such as mature T. distichum trees, were in 
many cases almost completely covered by the flood waters and subjected to greater water velocities, 
debris catching and greater stress (Figure 2, see Sigafoos 1964; Johnson et al. 2000). Flood induced 
changes in riparian forests are difficult to assess and recovery times are uncertain. It is known that high 
magnitude floods do recur at regular though potentially widely-spaced intervals in this ecosystem, and 
are apparently an important part of its dynamics. They may even be required if recruitment of early 
successional species is to occur in these communities. 

Recruitment of early successional species in these floodplain forests seems to require light gaps 
in the forest very much like the light gaps m other forests caused by senility, disease, insects, wind throw 
and fn*e (Brewer and Menitt 1978; Garwood et al. 1979; Putz et al. 1983; Begon et al. 2006; Van Auken 
and Bush 2013). High magnitude floods do in fact reduce cover and create gaps; they also recur at 
predictable, though widely-spaced, intervals in this ecosystem (Baker 1975; Kochel 1988), and are 
apparently an important force and part of the dynamics of this ecosystem. 
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Table 1. Flood induced changes in community mean basal area (m /ha) for 11 of the highest basal area 
riparian species found in the Guadalupe River riparian forest. Scientific names, common names, pre¬ 
flood, post-flood and % changes are presented. Total pre and post-flood mean basal areas are significantly 
different* (paired /-test, P<0.05). 


BASAL AREA (mVha) 

SPECIES 

NAME 

PRE-ELOOD 

POST-ELOOD 

% CHANGE 

Taxodium distichum 

Bald cypress 

61.1 

44.1 

-33 

Catya illinoinensis 

Pecan 

24.1 

22.0 

-8 

Ulmiis americana 

American elm 

3.6 

2.9 

-19 

Celtis laevigata 

Hackberry 

3.3 

2.4 

-27 

Platamis occidentalis 

Sycamore 

2.1 

1.9 

-10 

Morns rubra 

Red mulberry 

1.1 

0.3 

-63 

Acer negimdo 

Box-elder 

0.9 

0.8 

-12 

Ilex decidua 

Possum-haw 

0.6 

<0.1 

-86 

Cornus drummondii 

Dogwood 

0.1 

0.1 

0 

Ulnius crassifolia 

Cedar elm 

0.1 

0.1 

0 

Diospyros texana 

Persimmon 

<0.1 

<0.1 

0 

Others (20 species. 


1.6 

2.1 

+31 

see Table 3) 





Total* 


98.6±19.2 

76.5±18.3 

Mean = -22 


Table 2. Flood induced changes in community mean species density (plants/ha) for 11 of the highest 
density riparian species found. Pre-flood, post-flood and % changes are presented. Total pre- and post- 
flood mean densities are significantly different* (paired /-test, P<0.05). 


DENSITY (plants/ha) 

SPECIES 

PRE-ELOOD 

POST-ELOOD 

%CHANGE 

Carya illinoinensis 

229 

227 

-1 

Celtis laevigata 

141 

103 

-27 

Acer negundo 

139 

96 

-30 

Ulmus americana 

63 

46 

-27 

Taxodium distichum 

58 

45 

-22 

Cornus drummondii 

55 

46 

-26 

Ilex decidua 

43 

31 

-28 

Diospyros texana 

40 

29 

-27 

Platanus occidentalis 

13 

12 

-8 

Morus rubra 

15 

10 

-33 

Ulmus crassifolia 

12 

12 

0 

Others (20 species, see 

328 

176 

-46 

Table 3) 




Total* 

1124+171 

833+112 

Mean = -26 
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Table 3. Other speeies found in the Guadalupe River riparian forest ineluding eommon names and type of 
woody speeies. These speeies are low basal area and low density speeies. 


Species 

Common Name 

Vegetation Type 

Vitis m ustangensis 

Mustang grape 

Vine 

Sapindus driimmondii 

Soapberry 

Tree 

Parthenocissus quinquefolia 

Virginia creeper 

Vine 

P tele a trifoliata 

Skunk-bush 

Shrub 

Fraximis americana 

White ash 

Tree 

Cephalanthus occidentalis 

Buttonbush 

Vine 

Madura pom if era 

Osage orange 

Tree 

Juglans nigra 

Black walnut 

Tree 

Sideroxylon lanuginosum 

Gum bumellia 

Tree 

Prosopis glandulosa 

Honey mesquite 

Tree 

Melia azedarach 

Chinaberry 

Tree 

Toxicodendron rydbergii 

Western poison ivy 

Vine 

Quercus virginiana 

Live oak 

Tree 

Amorpha fruticosa 

Bastard indigo 

Shrub 

Callicarpa americana 

American beautyberry 

Shrub 

Juniperus asheii 

Ash juniper 

Tree 

Salix nigra 

Black willow 

Tree 

Ailanthus altissima 

Tree of heaven 

Tree 

Smilax hona-nox 

Cat-brier 

Vine 

Tilia ameracana 

Carolina basswood 

Tree 


Table 4. Flood indueed ehanges in mean species basal area (cm~/plant) for seven selected riparian species, 
including the four highest community total basal area species. Pre-flood, post-flood and % changes are 
presented. 


MEAN BASAL AREA 

SPECIES 

PRE-ELOOD 

POST-ELOOD 

% CHANGE 

Taxodium distichiim 

11,478 

9,450 

-18 

Caiya iUinomensis 

992 

845 

-15 

Ulmus americana 

878 

938 

+7 

Ulmiis crass if alia 

315 

271 

-14 

Celtis laevigata 

245 

247 

+1 

Corniis driimmondii 

13 

11 

-15 

Diospyros texana 

5 

6 

+20 
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Figure 1. Photographs of damage to the 
riparian forest of the Medina River. Top 
figure is ground level and shows the 
removal of the tops of all Taxodium trees 
on the river edge. The middle and bottom 
figures are aerial photographs showing a 
number of different kinds and amounts of 
damage to the floodplain forest. In most 
eases parts of the riparian forest were 
eompletely removed or destroyed by the 
flood. The middle pieture (aerial 
photograph), lower, eentral part shows all 
of the Taxodium trees with branehes 
removed and apparent total removal of any 
forest on the right side of the photograph. 
Bottom pieture (aerial photograph), lower 
eenter, next to the river shows Taxodium 
trees down or apparently washed away as 
the flood waters washed over the river 
edge. Photographs taken from Sullivan 
1983, with permission. 
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Figure 2. Examples of damage to Taxodium distichum and the riparian forest of the Guadalupe River after 
the 1978 Guadalupe River flood. Upper left, a mature Taxodium tree apparently undamaged with 
uprooted trees or broken stems on the bank of the river. Upper right, two mature Taxodium trees with tops 
removed. Lower left, piles of debris on the river bank with a large Taxodium tree standing in the 
approximate eenter with a large, broken limb eaught in the upper branehes approximately 16 m (63 ft.) 
above ground level. Lower right, large Taxodium tree with the two large top stems broken off and 
removed. Two broken or uprooted trees are at a diagonal in the photograph. Photographs were taken by 
the senior author. 
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Key to Juniperus occidentalis Hook, and J. o. forma corbetii. 

Robert P. Adams 

Biology Department, Baylor University, Box 97388, Waeo, TX 

76798, USA 

email Robert_Adams@baylor.edu 
and 

Mark Corbet 

7376 Southwest MeVey Ave., Redmond, OR 97756 

ABSTRACT 

Juniperus occidentalis forma corbetii R. P. Adams ean be readily distinguished by its shrubby 
habit and eompaet foliage vs. typieal J. occidentalis plants, that have a strong eentral axis and loose, open 
foliage (Fig. 2). A key is presented to aid in the identifieation of J. o. forma corbetii. Published on-line 
www.phytologia.orgP/zy/o/ogzdf 99(@).- 238-240 (Dec. 18, 2017). ISSN 030319430. 

KEY WORDS: Juniperus occidentalis, Juniperus occidentalis forma corbetii R. P. Adams, Key, 
identifieation, Cupressaeeae. 


Juniperus occidentalis, J. grandis (=J. occidentalis var. australis) and J. osteosperma are three 
very elosely related junipers in the western United States (Adams, 2014, Vasek, 1966). Reeently, a new 
taxon, J. o. forma corbetii R. P. Adams was deseribed (Adams, 2012) based on its shrubby form and 
eompaet foliage (Fig. 1,3) vs. typieal J. occidentalis plants, that have a strong eentral axis and loose, open 
foliage (Fig. 2). 



Figure 1. Mark Corbet with an extremely shrubby 
habit of J. occidentalis forma corbetii, 32 km east of 
Bend, OR (ef Adams 11949-11951) at the type 
loeality. 

Figure 2.(right). Juniperus occidentalis near Sage 
Hen, CA, with a strong eentral axis, and loose, open 
foliage. 
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Figure 3. Other plants of J. o. forma corbetii with less extreme shrubby habit. 

In addition, Adams (2012) reported the volatile leaf oil of J. occidentalis is dominated by 
sabinene, p-eymene, eitronellol and bomyl aeetate , whereas f corbetii shrubs, east of Bend, OR, have 
larger amounts of p-eymene (20.0) and bomyl aeetate (24.5%). 

Prineipal eoordinates analysis using 42 terpenoids showed (Adams, 2012) the f corbetii shmbs, 
east of Bend, were elearly separated (see Fig. 3, Adams, 2012). 

Overall, Adams (2012) found the leaf essential oils of populations of J. occidentalis were rather 
uniform exeept for the populations at the extremity of the range, and for the shmbby forma corbetii, east 
of Bend, OR. 

For identifieation of f corbetii from speeimens, it seems most effieient to utilize the foliage 
density and length of ultimate branehlets. This is illustrated in Fig. 4. Notiee that the foliage of J. 
occidentalis is open and individual branehlets are elearly visible (Fig. 4, left). In eontrast, the foliage of f 
corbetii is eompaeted, with short sub-branehing (Fig. 4, right). 



Figure 4. Comparison of the foliage of J. occidentalis and f corbetii. 

Key to J. occidentalis and f corbetii : 

1. Foliage loose and open, sub-branehlets elongated, trees with strong eentral axis. J. occidentalis 

1. Foliage eompaet and dense, sub-branehlets shortened, shmbs. J. o. forma corbetii 
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ABSTRACT 

A new speeies, Lycianthes textitlaniana E. Dean (Solanaeeae), is deseribed from the Sierra Madre 
del Sur of the state of Oaxaea, Mexieo. The affinities of L. textitlaniana are uneertain, beeause it differs 
from other known speeies of Lycianthes in its eombination of glandular pubeseenee, long ealyx 
appendages, and turbinate fruit with dry mesoearp. The seeds and pubeseenee of L. textitlaniana are 
somewhat similar to those of L. pringlei (B. E. Rob. & Greenm.) Bitter of western Mexieo. To help with 
identifieation of this new speeies, L. textitlaniana is eompared here with L. pringlei, L. surotatensis 
Gentry, L. gorgonea Bitter, L. purpusii (Brandegee) Bitter, and L. dejecta (Femald) Bitter, and a 
photograph of the seed is provided. Published on-line www.phvtologia.org Phytologia 99(4): 241-248 
(Dec. 18, 2017). ISSN 030319430. 

KEY WORDS: Lycianthes, Solanaeeae, Mexieo, Oaxaea, taxonomy 


Lycianthes (Dunal) Hassler (Solanaeeae) is native to both the New and Old Worlds and ineludes 
150 to 200 speeies (Hunziker, 2001). The majority of its taxa are distributed in the New World (from 
Mexieo to Argentina), with ea. 40 taxa native to Mexieo (ea. 15 endemie, one of them deseribed here) 
(Villasenor, 2016). The genus is the elosest relative of the ehile pepper genus Capsicum E. (Walsh & 
Hoot, 2001; Bohs & Olmstead, 1997). The German botanist Georg Bitter, who monographed the genus 
Lycianthes, was the first worker to point out that Lycianthes and Capsicum share a similar ealyx 
morphology in whieh the five sepal lobes are truneated into a sleeve, below whieh may protrude five to ten 
appendages (eommonly ealled ealyx teeth) (Bitter, 1919). However, while Capsicum speeies have anthers 
that dehisee by longitudinal slits, the speeies of Lycianthes typieally have porieidal anther dehiseenee. 

While working on speeies deseriptions for all Mexiean and Central Ameriean speeies of 
Lycianthes, the first author found a speeimen of a new Lycianthes speeies under the name Lycianthes 
gorgonea Bitter at the Herbario Naeional, Universidad Naeional Autonoma de Mexieo (MEXU). It was 
immediately obvious that this was a new speeies due to a unique eombination of eharaeters: glandular 
pubeseenee, long ealyx appendages, turbinate fruits with dry mesoearp, and seeds with deeply indented 
“pits” between wavy eell walls, produeing an ornamented surfaee (Figures 1 and 2). The turbinate fruits 
with dry mesoearp are unlike any other known Lycianthes. The pubeseenee and seed surfaee somewhat 
resemble those of L. pringlei, but it differs from that speeies in having very long ealyx appendages. 




242 


Phytologia (Dec 18, 2017) 99(4) 


This species is only known from one locality in Mexico and a single collection made in August, 
2006 by Alma Zarate-Marcos in a canyon by Colonia Nueva, to the west of the town of Santiago 
Textitlan, Sierra Madre del Sur, Oaxaca. Ms. Zarate-Marcos made the collection as part of a lloristic 
project managed by one of the co-authors (Sal as-Morales) in her role as Director of the “Sociedad para el 
Estudio de los Recursos Bioticos de Oaxaca” (SERBO). Santiago Textitlan is a community with forest 
management practices certified by the Forestry Stewai'dship Council and has been recognized for its good 
forest management (Blaclcman et af, 2015). Before this project, the flora of the region of Santiago 
Textitlan was undocumented. Since the initiation of the Santiago Textitlan floristic project, nearly 7500 
specimens in 140 families, 517 genera and 1172 species have been collected in the region (figures from 
the SERBO database). In addition, several new species fi'om this region have been described (Borhidi et 
af, 2008; Turner, 2013). 

Althougli we are describing this species based on very limited material (only the type collection), 
the material is very complete, with flowers, mature fruits, and seeds. By describing and naming this very 
unusual species of Lycianthes, we are adding to the knowledge of the biodiversity of Oaxaca, as well as 
facilitating the identification of future collections of this species. Termmology here follows that of 
previous papers on this genus (Dean, 1998; Dean et af, 2007). 

Lycianthes textitlaniana E. Dean, sp. nov. TYPE: Mexico: State of Oaxaca, Distrito Sola de Vega, 
Municipio Santiago Textitlan, Colonia Nueva, 1562 m, 18 Aug 2006, Alma Zarate Marcos AZM-274 
(holotype: MEXU-1229513 (Figure 1); isotype: SERO-45464). 

Diagnosis. This species is unlike any other Lycianthes species in its combination of glandular trichomes, 
long calyx appendages, dry turbinate and apiculate berries lacking fleshy mesocarp, and round-edged 
seeds with pitted surfaces. Indument and seeds are somewhat similar to those of L. pringlei, but L. 
pringlei lacks calyx appendages and has ovoid fruits with fleshy mesocarp. 

Shrub or vine, ca. 0.5 m tall (reported as an herb on the label of the only known collection, but stems are 
obviously woody). Indument of clear to white, multicellular, simple, glandular (glandular tip golden- 
yellow to grey), wavy to patent trichomes 0.25-1.5 mm long. Steins pale green (drying tan) when young, 
not much compressed when dried, becommg brown and woody with age, moderately to densely 
pubescent. Upper sympodial branching points a mixture of dichasial and monochasial branching (mostly 
monochasial), the upper sympodial units 1-5 cm long, 1-2.5 mm in diameter, difoliate, the leaves often 
geminate or unifoliate due to loss of small geminate leaves. Leaves simple, usually paired and unequal in 
size, the large ones with blades 3-7 x 1.5-3 cm, the smaller ones 1.5-5 x 0.5-2 cm, the leaf pairs similar 
in shape, the blades ovate, to elliptic, chartaceous, moderately to densely pubescent; major veins 4—5 on 
either side of tlie midvein; base tmncate, cuneate, or attenuate, sometimes oblique; margin irregularly 
wavy; apex acuminate; petioles absent or 0.1-1 cm long. Flowers 5-merous, solitary, axillary, the axes 
moderately to densely pubescent; peduncles absent; pedicels ca. 20-25 mm long and erect in flower, 20- 
30 mm long and arching to deflexed in fruit; calyx 3-3.5 imn long, 3.5-4 mm in diameter, campanulate, 
moderately to densely puberulent, the margin truncate, with 10 spreading, linear, basally flattened 
appendages, 7-10 mm long emerging ca. 0.5 mm below calyx rim; fruiting calyx accrescent, widely 
bowl-shaped, 4-6 mm long, 9-12 mm in diameter, the appendages to 15 mm long; corolla oriented 
horizontally, 1.2-1.5 cm long, rotate, mostly entire (with shallow notches), with abundant interpetalar 
tissue, purple adaxially, color of abaxial side unknown, glabrous, open when collected, uncertain of 
timing of diurnal movements; stamens unequal, with four short filaments 1-2 nun long and one long 
filament 2-3 mm long, glabrous, anthers 3.5-4 mm long, elliptic to lanceolate, not connivent, yellow, the 
pores ovate, nearly tennmal, dehiscing upward, not opening into longitudinal slits; ovary not exammed; 
style linear, ca. 7 mm long, glabrous; stigma capitate, decurrent down two sides. Fruit a dry berry lacking 
juicy mesocarp, 13-15 mm long, 9-12 mm in diameter, turbinate, the tip apiculate, pale greenish orange 
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(possible darker) when mature, glabrous, laeking selerotie granules; seeds ea. 50-80 per fruit, 1.75-2 x 
1.25-1.5 mm, round-edged, laterally compressed, reniform, tan-orange, the surface pitted (Figure 2). 

Distribution. Lycianthes textitlaniana is only known from the type collection from the state of 
Oaxaca, Mexico, in the District of Sola de Vega, just outside of the town of Santiago Textitlan at 1562 m 
in elevation (Figure 3). 

Habitat. Three of the co-authors (Dean, Torres-Colin, and Hinojosa-Espinosa) attempted to locate 
the new species in the field on September 7-8, 2017, after obtaining permission to explore the collection 
area from local authorities on September 7. Despite extensive exploration of the canyons near the 
latitude/longitude coordinates given on the collection label, the species was not located, although another 
Lycianthes species, L. acapulcensis (Baillon) D’Arcy, collected by Ms. Zarate-Marcos in the same 
location on the same date was encountered. The vegetation of the canyons that we explored consisted of 
Pine-Oak forest with elements of cloud forest and palms with the following species: Quercus 
peduncularis Nee, Arbutus xalapensis Kunth, Oreopanax, Senna hirsuta (L.) H. S. Imin & Bameby, 
Piper umbellatum L., Miconia glabrata Cogw., Psacalium cirsiifolium (Zucc.) H. Rob & Brettell, Bidens 
pilosa L., Podachaenium eminens (Lag.) Sch. Bip., Achimenes antirrhina (DC.) C. V. Morton, Crusea 
hispida (Mill.) Rob., Begonia gracilis Kunth, Hyptis capitata Jacq., Micropleura renifolia Lag., 
Phaseolus leptostachyus Benth,, and Canavalia hirsuta (M. Martens & Galeotti) Standi. The trunks of the 
trees in the areas had frre scars, and members of the community reported that a fire went through the area 
sometime after 2006. The effect of the fire on this species is not known but could be one reason that we 
could not locate it at the type locality. 

Phenology. The type specimens, collected in mid-August, each have one flower and several 
mature fhiits. Therefore, the species most likely begins flowering in July, near the beginning of the rainy 
season, and probably continues fruiting in September. 

Etymology. The speciflc epithet honors the community of Santiago Textitlan, Oaxaca. 

DISCUSSION 

Lycianthes textitlaniana is unlike any other known species in the genus, and it can be easily 
distinguished from other Mexican species with glandular trichomes {L. pringlei and L. surotatensis) by its 
turbinate fruit shape and dry mesocarp (Table 1). A number of Lycianthes species have calyx appendages 
as long as those found in L. textitlaniana, but none of the other species with long calyx appendages also 
have simple, glandular trichomes and turbinate fruit. Table 1 provides the characteristics of four other 
species with relatively long calyx appendages for comparison {L. surotatensis, L. gorgonea, and L. 
purpusii). Turbinate fhiits are found in several of the Mexican species of series Metonodontae (such as 
L. ciliolata (M. Martens & Galeotti) Bitter and L. dejecta (Femald) Bitter) (Bitter, 1919). The other 
species with turbinate fruits have eglandular trichomes. The seeds of L. textitlaniana are most similar to 
those of L. pringlei in having a pitted seed surface (Table 1; Dean et al., 2007). The closest relative of L. 
textitlaniana is not known. 
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Figure 1. Image of holotype of Lycianthes textitlaniana, Alma Zarate Marcos AZM-274 (MEXU- 
1229513). Use of speeimen eourtesy of the Herbario Naeional de Mexieo, Universidad Naeional 
Autonoma de Mexico. 
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Figure 2. Seed of Lycianthes textitlaniana (taken from the holotype Alma Zarate Marcos AZM-274 
(MEXU-1229513)); seale bar equals 0.3 mm. Photograph, by D. MeNair, produeed by staeking 45 images 
taken with a lOx objeetive lens attaehed to a Nikon DSLR. Use of speeimen eourtesy of the Herbario 
Naeional de Mexieo, Universidad Naeional Autonoma de Mexieo. 
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Figure 3. Map of known distribution of Lycianthes textitlaniana. 
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Table 1. Comparison of six characters in Lycianthes textitlaniana, L. pringlei, L. surotatensis, L. 
gorgonea, L. purpusii, and L. dejecta. 


Species 

textitlaniana 

pringlei 

surotatensis 

gorgonea 

purpusii 

dejecta 

Habit 

Probably 
shrub or vine 

Shrub 

Shrub 

Vine 

Vine 

Perennial 

herb 

Trichome 

type 

Simple, 

glandular 

Simple, 

glandular 

Simple, 

glandular 

and 

eglandular 

Simple, 

eglandular 

Simple to 
branched, 
eglandular 

Branched to 

stellate, 

eglandular 

Calyx 
appendage 
length in 
flower 

7-10 mm 
long 

Appendages 
reduced to 
protuberances 

2-10 mm 
long 

7-15 mm 
long 

8-12 mm 
long 

2-7 mm 
long 

Fruit 

shape 

Turbinate, tip 
apiculate 

Ovoid, tip 
rounded 

Spheric to 
depressed 
spheric, tip 
rounded 

Spheric to 
depressed 
spheric, tip 
rounded 

Spheric, tip 
rounded 

Turbinate, 
tip apiculate 

Mesocarp 

Scant pulp, 
appearing 
dry, seeds 
free 

Juicy pulp, 

surrounding 

seeds 

Juicy pulp, 

surrounding 

seeds 

Juicy pulp, 

surrounding 

seeds 

Juicy pulp, 

surrounding 

seeds 

Juicy pulp, 

surrounding 

seeds 

Seed 

surface 

Pitted 

Pitted 

Lacking pits 

Lacking pits 

Lacking pits 

Lacking 
pits, fibrils 
on walls 
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ABSTRACT 

Recently, Sobierajska et al. (2016), using nSSR and morphology, showed that Juniperus 
drupacea exliibited differentiation between Greece and Turkey/ Lebanon, suggestive of Pleistocene 
genetic isolation. Here, we report that leaf terpenoids and DNA sequence data support their hypothesis by 
confirming differentiation between Greece and Turkey/ Lebanon/ Israel. The leaf oils of the Turkey/ 
Lebanon plants contained one unique terpene (trans-verbenol, 0.1-1.4%) that was absent in the Greece 
plants. The Greece oil contained tliree terpenes not found in the Lebanon/ Turkey plants: (ar)-curcumene 
(2.2%), p-alaskene (0,3%) and a-alaskene (0.4%), Four other terpenes were in higher concentration in the 
Greece oils: camphene (0.4%), 5-3-carene (TO.9%), p-mentha-l,5-dien-8-ol, isomer (0.3%) and 4- 
terpineol (0.3%). Three terpenes were higher in Turkey and Lebanon oils: a-pinene (10.5 - 32.9%), 
hexadecanoic acid (0.4 - 1.4%) and trans-totarol (0.3 - 1.2%). Only one SNP was found (in iirDNA) that 
separated Greece from Turkey-Lebanon-Israel. No informative SNPs were found in petN-psbM, trnS- 
tmG, tmD-trnT or tniL-tmF cp regions. Published on-line www.phytologia.org Phytologia 99(2): 249- 
257 (Dec. 18, 2017). ISSN 030319430. 

KEY WORDS: Juniperus drupacea, Cupressaceae, geographic variation, terpenoids, nrDNA, cp DNA. 


Recently, Sobierajska et al. (2016) reported on variation in nSSR and moiphology fox Juniperus 
drupacea Labill. from Greece, Turkey and Lebanon. They found morphology clearly separated Greece 
plants from those in Turkey and Lebanon (Fig. 1). In addition, they found some differentiation between 
plants from the Taurus (Turkey) and Lebanon mountains (Fig. 1). 
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Analysis using nSSR data gave a similar 
pattern, with some differentiation within Greeee, 
and within Turkey populations. 


Sobierajska et al. (2016) discovered a 
lower level of diversity in nSSR data in the 
Greece populations versus more diversity in 
Turkey and Lebanon populations. They 
reasoned that genetic isolation during the 
Pleistocene had led to tire differentiation 
between the European and Asia populations. 


Figure 1. Clustering of J. drupacea populations by 
morphological data (adapted from Sobierajska et al. 
(2016)). 

Junipems drupacea is ancestral m the phytogeny of Juniperus (Adams, 2014, Adams and 
Schwarzbach, 2013). Adams (1997) reported the composition of the volatile leaf oil of J. drupacea was 
low in a-pinene in Greece (3.5-5.1%) but high in Turkey (14.3%) and the Crimea (cultivated) (22.1%) 
and high in limonene (+p-phellandrene) m Greece (46.7 - 48.4 %), Turkey (55.6%) and Crimea (50.3%). 
Juniperus drupacea oil was high, but variable m 5-3-carene (7.0 - 22.3%). The analysis by Adams (1997) 
based on a few samples from Greece and literature reports from Turkey and Crimea. 

The purpose of the present study was to present a more complete report on the composition of the 
volatile leaf oil of J. drupacea from Greece, Turkey and Lebanon, utilizing additional samples, all of 
which were analyzed in the same lab. In addition, nrDNA and cp DNA sequencing was perfonned to 
detennine if the differentiation reported by Sobierajska et al. (2016) is discemable in DNA sequencing 
(and in leaf volatile oils data). 



METHODS 

Specimens for DNA and volatile leaf oils were collected of J, drupacea from: 

Greece: 

2-3m tall, strong central axis, common with Quercus coccifera. Female cones falling off Polluiation 
in Sept? On highway, 18 Ion E of Tripolis at Achladokampos Pass, 700m, Greece, Prov. Arkadia. 7 
Aug. 1987, Coll: Robert P. Adams 5651,5652, 

small trees to 3 m, female common in the Achladokampos Pass but all trees young. 18 km E. of 
Tripolis. Prov. Arkadia. Eat. 37° 31.44'N, Long. 22° 37.92E., 560m, Greece, 29 Aug. 1999, Coll: 
Robert P. Adams 8795, 

large trees to 10m x 0.5 m dbh, male, abundant 7.3 km E of Kosmos. Pamonas Mts., Mt. Pamonas Ski 
iTft, 37° 08.00'N, 22° 45.63E, 870m, Greece, 29 Aug. 1999, Coll: Robert P. Adams 8796, 

pyramidal trees to 3-4m tall, strong central axis, common with Quercus coccifera, female cones 
falling off. Pollination in Sept? some trees monecious. On highway 7, 18km E of Tripoli at 
Achladokampos Pass, Country: Greece, Prov. Arkadia, 37° 31' 25" N, 22° 38' 19" E., 600m, 21 
May 2017, Coll. Robert P. Adams 15227-15236. 

Israel: 

Mt Hermon, 33.31° N, 35.77° E, 1700m, June 2015, Israel, Coll. Hagar V. Leschner sn. Lab Ace. 
Robert P. Adams 14580, 14581 

Lebanon: 
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Ehmej, on the road up the mountain. 34° 07' 58" N, 35° 49' 04" E, 1350 m, 22 Sept 2015, Eebanon, 
Coll. Madga Bou Dagher-Kharrat, Ehmej 1-5 (=EEB-2), Eab Aee: Robert P. Adams 14675-14679; 

Mrebbine, Sir Danniyeh, 34° 24' 29" N, 36° 05' 13" E, 1400 m, 20 Sept 2015, Eebanon, Coll. Madga 
Bou Dagher-Kharrat, Mrebbine 1-5 (=EEB-3, ~=Dufram), Eab Aee: Robert P. Adams 14680-14684. 

Mehati Keserwan, 34° 02' 02" N, 35° 45' 53" E, 1180 m, 22 Sept 2015, Eebanon, Coll. Madga Bou 
Dagher-Kharrat, Mehati Keserwan 1-5 (=EEB-1), Eab Aee: Robert P. Adams 14685-14689 
Turkey: 

Camliyayla (Namrun), above Kole village, vieinity of Hotel Dag, 37° 10' 6.41" N, 34° 36' 3.8" W., 
1250 m, Jan. 2017; Turkey, Coll. S. Tugrul Koruklu sn, Eab Aee. Robert P. Adams 15220; 

Alanya, Derekoy, Gedevet Pasture, Park Orman. 36° 39' N, 32° 10' E, elev. ea 1200 m, 19 Sept 2015, 
Prov. Antalya: Turkey, Coll. Tugrul Mataraci & Sezai Vzun 2015-21, Eab Aee: Robert P. Adams 
14693; 

Gedevet Oastyrem Derinee, on road to Sapaeietlik, in the river bed. 36° 38' 23.27" N, 32° 03' 38.44" 
E., 1036 m, 13 Sep. 2017, Prov. Antalya, Turkey, Coll. T. Mataraci, O. Demir & B. Cingay 8332, 
Eab Aee. Robert P. Adams 15307 

Gedevet Oastyrem Derinee, on road between Antalya and Askeki eounty,30 km, 37° 03' 32.85" N, 31° 
44' 47.34" E, 1132m, 12 Sep. 2017, Prov. Antalya, Turkey, Coll. T. Mataraci, O. Demir & B. Cingay 
BCNG 8328A, Eab Aee. Robert P. Adams ,15308 

Voueher speeimens are deposited in the herbarium, Baylor University (BAYEU). 

One gram (fresh weight) of the foliage was plaeed in 20 g of aetivated siliea gel and transported 
to the lab, thenee stored at -20T C until the DNA was extraeted. DNA was extraeted from juniper leaves 
by use of a Qiagen mini-plant kit (Qiagen, Valeneia, CA) as per manufaeturer's instruetions. 

Amplifieations were performed in 30 pi reaetions using 6 ng of genomie DNA, 1.5 units Epi¬ 
centre Fail-Safe Taq polymerase, 15 pi 2x buffer E (petN-psbM), D (maldehy) or K (nrDNA) (final 
eoneentration: 50 mM KCl, 50 mM Tris-HCl (pH 8.3), 200 pM eaeh dNTP, plus Epi-Centre proprietary 
enhaneers with 1.5 - 3.5 mM MgCC aeeording to the buffer used) 1.8 pM eaeh primer. See Adams, Bartel 
and Priee (2009) for the ITS and petN-psbM primers utilized. The primers for tmD-tmT, tmE-tmF and 
tmS-tmG regions have been previously reported (Adams and Kauffmann, 2010). The PCR reaetion was 
subjeeted to purifieation by agarose gel eleetrophoresis. In eaeh ease, the band was exeised and purified 
using a Qiagen QIAquiek gel extraetion kit (Qiagen, Valeneia, CA). The gel purified DNA band with the 
appropriate sequeneing primer was sent to MeEab Ine. (San Franeiseo) for sequeneing. Sequenees for 
both strands were edited and a eonsensus sequenee was produeed using Chromas, version 2.31 
(Teehnelysium Pty Etd.) or Sequeneher v. 5 (geneeodes.eom). Sequenee datasets were analyzed using 
Geneious v. R7 (Biomatters, available from http://www.geneious.eom/) , the MAFFT alignment program. 
Further analyses utilized the Bayesian analysis software Mr. Bayes v.3.1 (Ronquist and Huelsenbeek 
2003). For phylogenetie analyses, appropriate nueleotide substitution models were seleeted using 
Modeltest v3.7 (Posada and Crandall 1998) and Akaike's information eriterion. 

Fresh leaves (200 g) were steam distilled for 2 h using a eireulatory Clevenger-type apparatus 
(Adams, 1991). The oil samples were eoneentrated (ether trap removed) with nitrogen and the samples 
stored at -20°C until analyzed. The extraeted leaves were oven dried (100°C, 48 h) for determination of 
oil yields. 

The oils were analyzed on a HP5971 MSD mass speetrometer, sean time 1/ see., direetly eoupled 
to a HP 5890 gas ehromatograph, using a J & W DB-5, 0.26 mm x 30 m, 0.25 mieron eoating thiekness, 
fused siliea eapillary eolumn (see Adams, 2007 for operating details). Identifieations were made by 
library searehes of the Adams volatile oil library (Adams, 2007), using the HP Chemstation library seareh 
routines, eoupled with retention time data of authentie referenee eompounds. Quantitation was by FID on 



252 


Phytologia (Dec. 18, 2017) 99(4) 


an HP 5890 gas chromatograph using a J & W DB-5, 0.26 mm x 30 m, 0.25 micron coating thickness, 
fused siliea eapillary eolumn using the HP Chemstation software. 

RESULTS 

The eomposition of the leaf essential oils from Greeee (avg. 10 plants), Turkey (3 plants), and 
three loeations in Lebanon: Leb-Mb (avg. 5 plants), Leb-Me (avg. 5 plants), and Leb-Eh (avg. 5 plants) 
are given m Table 1. The leaf oils of the Turkey/ Lebanon plants contamed one unique terpene (trans- 
verbenol, 0.1-1.4%) that was absent in the Greeee plants (Table 1). The Greeee oil eontained three 
terpenes not found in the Lebanon/ Turkey plants; (ar)-curcumene (2.2%), p-alaskene (0.3%) and a- 
alaskene (0.4%) (Table 1). Four other 
terpenes were in higher eoncentration in the 
Greece oils: camphene (0.4%), 8-3-earene 
(10.9%), p-mentha-l,5-dien-8-ol, isomer 
(0.3%) and 4-terpineol (0.3%). Three 
terpenes were higher (Table 1) in Turkey and 
Lebanon oils: a-pinene (10.5 - 32.9%), 
hexadeeanoie aeid (0.4 - 1.4%) and trans- 
totarol (0.3 - 1.2%). 

Factoring the similarity matrix (based 
on quantitative comparisons of sixteen major 
eompounds), resulting in thi*ee eigenroots that 
aeeounted for 84% of the variance among 
OTUs. The eigenroots appeared to asymptote 
after the third eigenroot. Principle Coordinate 
Ordination (PCO) analysis of the similarity 
matrix revealed the major trend (eoordinate 1, 

40%, Fig. 2) was the differentiation of the 
Greece oil from the oils of Turkey and 
Lebanon (Fig. 2). 

The seeond eoordinate (28%) separated 
Lebanon (Elimej, Mchati, Mrebbine) and T1 
(Turkey, 14693) from two other Turkey 
individuals (T2, T3, Fig. 2). 

The third eoordinate (17%) separated T2 from T3 and (Eh, Me) from (Mb, Tl, Fig. 2). 

Plotting eontoured similarity elustering on a geographic map gives one the spatial aspect of the 
variation (Fig. 3). Notice that all the Turkey-Lebanon oils cluster at a 0.75 similarity (Fig. 3), whereas the 
Greece oil is only admitted into the final eluster at a similarity of 0.52. This shows the magnitude 
(eonsiderable) of differentiation between Greece and Turkey-Lebanon oils. In general, the clustering is 
similar to that of the morphology (ef Fig. 1 vs. Fig. 3). 

It is also interesting that two of the Lebanon populations (Eh, Me) eluster at 0.93, as shown by 
their nearly identieal pereent eoneentrations of eomponents in Table 1. One Eebanon population (Mb) 
elusters with a Turkey plant (Tl) at 0.92 (also notiee the similarity in their oils in Table 1). 



Figure 2. PCO of Greeee population (avg. 10), three 
Eebanon populations (5 each, Ehmej, Mchati, 
Mrebbme) and three Turkey plants (Tl, T2, T3). 
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Figure 3. Contoured elustering based 
on 16 major terpenoids. See text for 
diseussion. 


It is informative to eompare the morphologieal 
elustering (Sobierajska et al. 2016) versus clustering 
using 16 terpenoids (Fig. 4). Both morphology and 
terpenoids show major differentiation between Greece 
(Europe) and Asia plants. Both data show an individual 
from Lebanon clustering with a plant from Turkey 
(Morphology: Leb-2, TAT-1; Terpenoids: Tl, LEB3- 
Mrebbine, Fig. 4). 


Among the Turkey plants, morphology displays 
less diversity than is seen in terpenoids (Fig. 4). 
Otherwise, the clustering is very similar between the 
morphology and terpenoids data sets. 


Figure 4. Comparison of clustering based on 
morphology (Sobierajska et al. 2016) and 
terpenoids (this study). 

In addition to the analysis of terpenoids, we also sequenced nrDNA (ITS) and four cp regions: 
petN-psbM, tmS-trnG, triiD-tniT and tniL-tniF, These gene regions produced 3,423 bp of data. Only 
one SNP was found (in mDNA) that separated Greece from Turkey-Lebanon-lsrael. No informative 
SNPs were found in petN-psbM, tmS-tmG, tmD-tmT or tniL-tniF cp regions. 

A Bayesian tree shows (Fig. 5) J. dnipacea (monotypic in sect. Caryocedrus) well resolved from 
Juniperus sect Jimiperus as previously shown (Adams, 2014, Adams and Schwarzbach, 2013). All the 
Greece plants group together and that clade is separated by one SNP in the nrDNA region. One SNP (in 
nrDNA) resolved the clade contauiing: Israel, 14580; Me Leb, 14685 and Me Leb, 14686. Eh Leb, 14686 
contained a unique SNP, but otherwise there was no variation found. Although the divergence of the 
Greece plants was based on only one SNP, it was consistent among the seven plants and supports the 
morphology, nSSR, and terpenoids in showing differentiation between Greece (Europe) and Asia plants. 
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Bayesian tree, 3,423 bp 

nrDNA, petN-psbM, trnS-G, trnD-T, trnL-F 

- 1 




Hesperocyparis arizooica 


Hesperocyparis bakeri 

Turkey. 15220 


Juniperus drupacea 
Beet. Caryocedrus 


- 7 ^ 


0.99. 


0.99. 




Turkey, 14693 
Israel. 14581 
Mb Leb, 14680 
Mb Leb, 14681 
- Eh Leb. 14676 
-Eh Leb, 14681 
Israel, 14580 
Me Leb, 14685 
Me Leb. 14686 
Greece. 5651 
Greece. 5652 

4 Greece. 8795 
Greece, 8796 
Greece, 15227 
Greece. 15228 
Greece. 15229 


Juniperus sect. Juniperus 
14 species 


Figure 5. Bayesian tree based on 
3,423 bp, nrDNA, petN-psbM, tmS- 
tmG, tniD-tmT and tmL-trnF. 
Numbers at the branch points are 
posterior probabilities. 


This study reconfirms the work of Sobierajska et al. (2016) that differentiation has occurred 
between Europe (Greece) and Asia (Turkey and Lebanon) in J. drupacea. Although the differences do 
not support the recognition of a new variety or subspecies, they present significant evidence that supports 
the thesis that the Aegean Sea presented a significant barrier during the Pleistocene to gene flow between 
Europe and Asia populations of J. drupacea and has contributed to the differentiation documented in this 
study and that of Sobierajska et al. (2016). 

A hypothetical ancestor of the species, Juniperus drupacea-pliocenica Rer., probably diverged 
prior to 30-35 Mya (Mao et al., 2010) and evolved on the plate of the European continent (Palamarev, 
1989; Kvacek, 2002; Palamarev, Bozukov, Uzunova, Petkova, Kitanov 2005). The colonization of 
terrains of the contemporaiy^ Peloponnese, Anatolia and the north-western part of the Arabian plateau 
might have occurred during Oligocene and the early Miocene (Burdigalian, 20.44—15.97 Mya), when a 
solid land comiection supporting plant migrations was formed between these regions (Rogl, 1999; 
Meulenkamp & Sissingh, 2003; Popov et al., 2006). Subsequently, the separation of the Pelopomiese 
microplate from Anatolia (Tortonian 11.60-27.24 Mya) created the Aegean Sea. This presented a barrier, 
triggering their differentiation, afterward magnified during Messinian salinity crisis. The geographic 
range of J. drupacea was reduced during Pliocene climate cooling and drying (Thompson, 2005). These 
climate changes probably caused altitudinal migration, resulting in spatial isolation betu^een populations 
inhabiting particular mountain massifs and, consequently, in subsequent genetic divergence (Sobierajska 
et al. 2016). The very dynamic clunate fluctuations durmg Pleistocene, with cold (glacial) and hot 
(interglacial) periods, additionally modified by humid and/or dty phases (Roberts et al. 2011; Dean et al. 
2015) influenced geographic distribution of organisms, including J. drupacea. These factors seem likely 
responsible for the current geographic range and at least partly for genetic, biochemical and 
morphological differentiation in populations of J. drupacea. 

ACKNOWLEDGEMENTS 

This research was supported in part with funds from Baylor University. Thanks to Amy TeBeest 
for lab assistance. 




Phytologia (Dec. 18, 2017) 99(4) 


255 


LITERTURE CITED 

Adams, R. P. 1991. Cedarwood oil - Analysis and properties, pp. 159-173. in: Modem Methods of Plant 
Analysis, New Series: Oil and Waxes. H.-F. Linskens and J. F. Jackson, eds. Springier- Verlag, 
Berlin. 

Adams, R. P. 2007. Identification of essential oil components by gas chromatography/ mass 
spectrometry. 4th ed. Allured Publ., Carol Stream, IL. 

Adams, R. P. 2014. The junipers of the world: The genus Jiiniperus. 4th ed. Trafford Publ., Victoria, BC. 

Adams, R. P. 1997. Comparisons of the volatile leaf oils of Jumperus drupacea Labill. from Greece, 
Turkey and the Crimea. JEOR 9:541-544. 

Adams, R. P. J. A. Bartel and R. A. Price. 2009. A new genus, Hesperocyparis, for the cypresses of the 
new world. Phytologia 91: 160-185. 

Adams, R. P. and M. E. Kauffmaim. 2010. Geographic variation in the leaf essential oils of Juniperus 
and comparison with J. occidentalis mdJ. osteosperma. Phytologia 92: 167-185. 

Adams, R. P. and A. E. Schwarzbach. 2013. Phylogeny of Juniperus usmg nrDNA and four cpDNA 
regions. Phytologia 95: 179-187. 

Dean JR, Jones MD, Leng MJ, Noble SR, Metcalfe SE, Sloane HJ, Sahy D, Eastwood WJ. 2015. Eastern 
Mediterranean hydroclimate over the late glacial and Holocene. Reconstmcted from the sediments of 
Nar lake, central Turkey, using stable isotopes and carbonate mineralogy. Quaternary Science 
Reviews 124: 162-174. 

Kvacek Z. 2002. A new juniper from the Paleogene of Central Europe. Feddes Repeitorium 113: 492- 
502. 

Mao K, Hao G, Liu J, Adams RP, Milne RI. 2010. Diversification and biogeography of Junipems 
(Cupressaceae): variable diversification rates and multiple intercontinental dispersals. New 
Phytologist 188: 254—272. 

Meulenkamp JE, Sissingh W. 2003. Tertiaiy palaeogeography and tectonostratigraphic evolution of the 
northern and southern Peri-Tethys platforms and the intennediate domains of the African-Eurasian 
convergent plate boundary zone. Palaeogeography, Palaeoclimatology, Palaeoecology 196: 209-228. 

Palamarev E. 1989. Paleobotanical evidences of the Tertiary histoiy and origin of the Mediterranean 
sclerophyll dendroflora. Plant Systematics and Evolution 162: 93-107. 

Palamarev E., Bozukov V., Uzunova K., Petkova A., Kitanov G. 2005. Catalogue of the Cenozoic plants 
of Bulgaria (Eocene to Pliocene). Phytologia Balkanica 11 (3): 215-364. 

Popov SV, Shcherba IG, Ilyina LB, Nevesskaya LA, Paramonova NP, Khondkarian SO, Magyar 1. 2006. 
Late Miocene to Pliocene palaeogeography of the Paratethys and its relation to the MediteiTanean. 
Palaeogeography, Palaeoclimatology, Palaeoecology 238: 91-106. 

Posada, D. and K. A. Crandall. 1998. MODELTEST: testing the model of DNA substitution. 
Bioinfonnatics 14: 817-818. 

Ronquist, F. and J. P. Huelsenbeck. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed 
models. Bioinformatics 19: 1572-1574. 

Sobierajska, K., K. Boratyhska, A. Jasihska, M. Dering, T. Ok, B. Douaihy, M. B. Daugher-Kharrat, A. 
Romo and A. Boratyhski. 2016. Effect of the Aegean Sea barrier between Europe and Asia on 
differentiation in Juniperus drupacea (Cupressaceae). Bot. J. Linnean Soc. : 1 -21. 

Roberts N, Eastwood WJ, Kuzucuoglu C, Fiorentino G, Caracuta V. 2011. Climatic, vegetation and 
cultural change in the eastern Mediterranean during the mid-Holocene enviromnental transitions. 
Holocene (Special Issue) 21: 147-162. 

Rogl F. 1999. Mediterranean and Paratethys. Facts and hypothesis of an Oligocene to Miocene 
paleogeography (short overview). Geologica Caipatica 50: 339-349 

Thompson JD. 2005. Plant evolution in the Mediterranean. Oxford: Oxford University Press. 



256 


Phytologia (Dec. 18, 2017) 99(4) 


Table 1. Leaf essential oil eompositions for Juniperus drupacea from Greeee, Turkey and Lebanon. 


921 

Tricyclene 

t 

t 

t 

t 

t 

t 

t 

924 

a-thujene 

t 

t 

t 

t 

t 

t 

t 

932 

a-pinene 

5.0 

10.5 

16.6 

12.5 

12.3 

32.9 

32.8 

945 

a-fenchene 

t 

t 

t 

0.1 

t 

t 

t 

946 

Camphene 

0.4 

t 

t 

0.1 

t 

0.2 

0.2 

969 

Sabinene 

0.2 

t 

t 

0.2 

0.2 

0.2 

0.2 

974 

P-pinene 

0.3 

0.9 

0.7 

1.0 

0.8 

1.8 

1.5 

988 

Myrcene 

2.3 

3.2 

2.7 

3.1 

2.6 

3.2 

3.0 

1002 

a-phellandrene 

t 

0.2 

t 

0.1 

t 

t 

t 

1008 

6-3-car ene 

10.9 

- 

1.2 

0.7 

0.4 

t 

t 

1014 

a-terpinene 

t 

t 

t 

t 

t 

t 

t 

1020 

p-cymene 

0.2 

t 

0.5 

0.1 

t 

t 

t 

1024 

Limonene 

31.3 

26.2 

17.4 

34.6 

32.0 

21.5 

21.3 

1025 

p-phellandrene 

28.1 

26.2 

17.4 

34.7 

31.2 

21.5 

21.3 

1054 

y-terpinene 

0.1 

t 

t 

0.1 

t 

t 

t 

1086 

Terpinolene 

1.1 

0.3 

0.6 

0.4 

0.4 

0.7 

0.5 

1099 

Linalool 

t 

0.2 

0.8 

0.1 

0.1 

0.2 

0.2 

1100 

n-nonanone 

- 

t 

t 

t 

t 

t 

t 

1118 

cis-p-iTienlha-2-en-1 -ol 

t 

t 

t 

t 

t 

- 

- 

1123 

a-camphenal 

t 

0.1 

0.5 

- 

0.2 

0.4 

0.5 

1132 

cis-limonene o.\ide 

0.1 

- 

- 

- 

- 

- 

- 

1135 

trans-pinocarveol 

- 

t 

0.9 

- 

0.2 

0.6 

0.7 

1136 

trans-p-mentha-2-en-1 -ol 

- 

0.1 

t 

0.1 

- 

- 

- 

1139 

trans-limonene oxide 

- 

- 

- 

- 

- 

- 

- 

1140 

trans-verbenol 

- 

0.2 

1.4 

0.1 

0.4 

0.7 

0.8 

1141 

Camphor 

t 

- 

- 

- 

- 

- 

- 

1156 

p-mentha-I,5-dien-8-ol, isomer 

0.3 

t 

t 

t 

t 

0.1 

t 

1166 

p-mentha-1,5-d ien-8-ol 

0.2 

t 

t 

t 

0.1 

0.2 

0.2 

1174 

terpmeii-4-ol 

0.3 

0.1 

t 

0.1 

0.1 

0.1 

0.2 

1176 

m-cymen-8-ol 

0.1 

t 

t 

t 

t 

t 

t 

1179 

p-cymen-8-ol 

0.1 

t 

t 

t 

t 

t 

t 

1186 

a-tcrpineol 

0.4 

0.3 

0.5 

0.2 

0.1 

0.1 

0.1 

1191 

cis-dihydrocarvone 

t 

- 

- 

- 

- 

- 

- 

1195 

Myrtenol 

t 

t 

0.2 

t 

t 

0.1 

0.2 

1195 

cis-piperitol 

t 

- 

- 

- 

- 

- 

- 

1204 

Irans-dihydrocarvone 

- 

- 

- 

- 

- 

- 

- 

1204 

Verbenone 

0.1 

t 

0.6 

t 

0.1 

0.2 

0.2 

1215 

trans-carveol 

0.1 

t 

0.3 

t 

0.3 

0.2 

0.3 

1223 

Cilronellol 

t 

- 

- 

- 

- 

- 

- 

1226 

cis-carveol 

t 

- 

- 

- 

- 

- 

- 

1239 

Carvone 

0.1 

0.2 

0.3 

0.1 

0.1 

t 

t 

1241 

carvacrol, methyl ether 

- 

t 

t 

t 

0.1 

t 

t 

1244 

car-3-en-2-one 

t 

- 

- 

- 

- 

- 

- 

1249 

Piperitone 

t 

- 

- 

- 

- 

- 

- 

1287 

bornyl acetate 

t 

t 

t 

t 

t 

t 

t 

1345 

a-cubebene 

0.5 

t 

0.4 

0.1 

0.2 

0.2 

0.2 

1374 

a-copaene 

0.5 

t 

0.7 

0.1 

0.4 

0.4 

0.5 

1387 

P-bourbonene 

t 

t 

t 

t 

t 

t 

t 

1417 

(E)-caryophyllene 

0.4 

0.6 

1.2 

02 

0.9 

0.9 

0.9 

1429 

cis-thujopsene 

t 

- 

- 

- 

- 

- 

- 

1448 

cis-muurola-3,5-diene 

t 

t 

t 

t 

- 

- 

- 

1452 

a-humulene 

t 

0.5 

0.7 

0.2 

0.4 

0.4 

0.4 

1465 

cis-muurola-4{ 14),5-diene 

0.3 

0.5 

0.5 

0.2 

t 

t 

t 

1478 

y-muurolene 

t 

0.7 

1.4 

0.1 

0.3 

0.3 

0.3 

1480 

gemiacrene D 

1.5 

2.9 

3.5 

1.5 

1.9 

1.9 

2.1 

1480 

(ar-)curciimene 

2.2 

- 

- 

- 

- 

- 

- 

1496 

Valencene 

- 

0.3 

0.4 

t 

0.1 

t 

0.1 

1498 

P-alaskcnc 

0.3 

- 

- 

- 

- 

- 

- 

1500 

a-muurolene 

0.2 

t 

0.4 

t 

0.2 

0.2 

0.2 

1505 

(E,E)-a-farnesene 

- 

- 

- 

- 

0.2 

t 

t 

1512 

a-alaskene 

0.4 

- 

- 

- 

- 

- 

- 

1513 

y-cadinene 

0.6 

1.3 

1.8 

t 

1.1 

0.6 

0.8 

1521 

trans-calamenene 

0.4 

t 

0.7 

0.1 

0.1 

0.1 

0.2 

1522 

5-cadinene 

t 

t 

0.7 

0.1 

0.1 

0.2 

0.3 
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KI 

Compound 

Greece 

avg. 

rurk-T2 

15307 

rurk-T3 

15308 

Turk-Tl 

14693 

Leb-Mb 

14691 

Leb-Mc 

14692 

Leb-Eh 

14690 

1533 

10-epi-cubenol 

0.2 

0.4 

0.5 

0.1 

0.2 

0.1 

0.1 

1533 

trans-cadina-1,4-diene 

0.3 

0.4 

0.5 

0.1 

0.2 

0.1 

0.2 

1561 

(E)-nerolidol 

- 

0.1 

0.4 

0.2 

0.1 

- 

t 

1570 

Dendrolasin 

- 

0.1 

0.1 

t 

- 

- 

- 

1582 

(ar-) tumerol 

t 

- 

- 

- 

- 

- 

- 

1582 

car^mphyllene o.x.ide 

0.3 

1.0 

1.5 

t 

0.7 

0.6 

0.6 

1586 

gennacra-4( 15),5,10( 14)- 
trien-l-al, isomer I 

t 

0.5 

0.5 

t 

0.3 

0.2 

0.2 

1594 

salvial-4( 14)-en- l-one 

t 

0.5 

0.7 

t 

0.3 

0.2 

0.2 

1608 

luimulene epo.xide 11 

0.4 

0.5 

0.8 

t 

0.3 

0.2 

0.3 

1610 

saIvial-4(14)-cn-l-one, isomer 1 

0.6 

0.9 

1.2 

0.1 

0.3 

0.2 

0.2 

1638 

epi-rx-cadinol 

0.2 

0.3 

0.4 

0.1 

0.3 

0.2 

0.3 

1640 

epi-a-muurolol 

0.2 

0.3 

0.3 

0.1 

0.3 

0.2 

0.3 

1652 

a-cadinol 

1.1 

1.8 

2.1 

0.7 

0.8 

0.6 

0.7 

1685 

germacra-4( 15),5,10( 14)- 
trien-l-al 

0.3 

0.7 

0.5 

0.2 

0.4 

0.3 

0.3 

1687 

eudesma-4( 15),7-dien-l-P-ol 

0.4 

0.7 

0.6 

0.3 

0.4 

0.3 

0.3 

1711 

Pentadecanal 

t 

t 

0.4 

t 

t 

t 

t 

1840 

hexahydrofarnesyl acetone 

- 

t 

t 

t 

t 

t 

t 

1959 

hcxadccaaoic acid 

t 

1.4 

1.3 

0.7 

0.4 

0.5 

0.4 

1987 

manoyl oxide 

1.3 

2.5 

4.5 

0.3 

4.0 

1.9 

0.8 

2055 

Abietatriene 

t 

0.5 

0.5 

0.2 

0.1 

0.3 

t 

2105 

Isoabicnol 

0.3 

1.3 

2.1 

0.6 

0.7 

0.8 

0.5 

2112 

phytol, isomer 1 

0.2 

1.3 

0.8 

0.4 

0.4 

0.6 

0.3 

2231 

methy l communate, isomer I 

- 

t 

t 

t 

0.1 

0.3 

0.2 

2314 

trans-totarol 

t 

0.3 

1.2 

0.3 

0.3 

0.5 

0.3 

2331 

trans-ferruginol 

t 

0.5 

0.6 

0.2 

t 

0.2 

t 


KI = linear Kovats Index on DB-5 column. Compositional values less than 0.1% are 
denoted as traces (t). Unidentified components less than 0.5% are not reported. 




